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Dear Sir: 

1. This declaration is to establish the status of the invention in the above- 
captioned U.S. Patent Application prior to March 26, 1993, the filing date of U.S. Patent 
No. 5,368,026, issued on November 29, 1994, to Swedlow et al. ("Swedlow"). 

2. We understand that Claims 1-10, 12-16, 18-25, 27-36, and 38-46 of the 
above-captioned U.S. Patent Application Ser. No. 10/677,050 (the '050 Patent 
Application), which was filed September 30, 2003, titled "Signal Processing Apparatus," 
have been rejected in the February 27, 2008 Office Action as being unpatentable over 
Swedlow. 

3. We are inventors of the subject matter of the above-referenced claims in 
the '050 Patent Application. 

4. At least prior to March 26, 1993, we conceived of and reduced to practice 
a signal processing apparatus such as the device and method described and claimed in 
the above-referenced claims of the '050 patent application. 

5. As evidence thereof, we provide herewith Exhibit A which is a design 
note ("the Design Note") describing the subject matter of the above-referenced claims. 
The Design Note was created and is dated prior to March 26, 1993. We conceived the 



ideas referenced therein prior to the creation of the Design Note. Any references to 
dates on the Design Note have been redacted. 

6. At Masimo Corporation ("Masimo"), the term "Stork" refers generally to 
Masimo's pulse oximetry technologies. 

7. At least prior to March 26, 1993, we compiled, loaded onto hardware, and 
used a software application called "STORK" for its intended purpose, that is, to measure 
blood oxygen saturation. 

8. As evidence thereof, we provide herewith Exhibits B and C which are 
screen captures of a directory of files, including B .c" source code files, ".h" header code 
files, and ".gdd" debugging code files compiled into the application file STORK. As 
shown on the screen captures, each of these files includes a last edited date prior to 
March 26, 1993, which has been redacted herein. 

9. Exhibits D and E are the header and source code files, titled "vector.h," 
and "vector.c," respectively, written in the C programming language, showing the 
mathematical functions for calculating a vector cross correlation. 

10. Exhibits F and G are the header and debugging files, titled "stork.h," and 
u stork.gdd," respectively, written in the C programming language. The header file 
defines the variable for cross correlation and the debugging file displays the value of 
cross correlation for instructive purposes. 

11. Exhibit H is the source code file, titled "iawa.c," written in the C 
programming language, which calculated the cross correlations. 

12. Exhibit I is the source code file, titled "satfil.c" written in the C 
programming language, which qualifies bins of signal data responsive to intensity 
signals from a detector monitoring a patient, the qualification including determining 
whether the calculated cross correlation is within expected limits. 

13. Exhibit J is the source code file, titled "ui.c," written in the C programming 
language, which is the user interface for the STORK application and which indicates an 
expected limit of cross correlation. 

14. The foregoing files and other exhibits demonstrate that we constructed 
and tested an embodiment and performed processes meeting every element of the 
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above-referenced claims prior to March 26, 1993 sufficient to demonstrate that the 
constructed embodiment operated for its intended purpose. 

15. We hereby declare that acts leading to the reduction to practice of the 
subject matter claimed in the '050 Patent Application were performed in the United 
States. 

i declare that all statements made herein of my own knowledge are true and that 
all statements made on information and belief are believed to be true, and further that 
these statements were made with the knowledge that willful, false statements and the 
like so made are punishable by fine or imprisonment, or both, under Section 1001 of 
Title 18 of the United States Code and that such willful, false statements may jeopardize 
the validity of the application or any patent resulting therefrom. 



Dated: By: 

Mohamed K. Diab 



I declare that all statements made herein of my own knowledge are true and that 
all statements made on information and belief are believed to be true, and further that 
these statements were made with the knowledge that willful, false statements and the 
like so made are punishable by fine or imprisonment, or both, under Section 1001 of 
Title 18 of the United States Code and that such willful, false statements may jeopardize 
the validity of the application or any patent resulting therefrom. 



Dated: Au^js-tJZ, Z££S_ By: A>J^L£&U 7z2 ^Z^^. 

Walter M. Weber 
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I declare that all statements made herein of my own knowledge are true and that 
all statements made on information and belief are believed to be true, and further that 
these statements were made with the knowledge that willful, false statements and the 
like so made are punishable by fine or imprisonment, or both, under Section 1001 of 
Title 18 of the United States Code and that such willful, false statements may jeopardize 
the validity of the application or any patent resulting therefrom. 



Dated: By: 

Massi E. Kiani 
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MIIK10 CORRELATION DETECTION n^i 9 " 
™ IU ™ ALGORITHM nDtes 


Design Note 
37 


Revision 

A 


File 

SDN37A.DOC 


Stork 


Page 

lof2 


Walter Weber 


Revision Date 


Original Dale 







The goal of this note is to propose an artifact detector for the Stork algorithm. Specifically, the artifact 
detector is designed such that it assigns a real number to each snapshot of data that it processes, and then 
compares that real number to a predetermined threshold value. Snapshots of data which yield artifact 
detector values that lie below the threshold value are said to contain significant amounts of artifact. 
Snapshots of data that yield artifact detector values which lie above the threshold value are said to contain 
insignificant amounts of artifact. 



In principle, the Stork algorithm accepts as input f 
bandpass) 



> signals which are of the form (after log and 



S a (0 = s a (t) + n a (t) 
S b (0 = s b (t) + n b (t) 



s b (t) = o>, s a (t) 

(2) 

n b(0 = ffl2naO- 

Any input signals which are not expressable in the above form are said to consist of artifacts. 
Consequently to determine whether or not a signal is arlifact free, it is sufficient to determine whether or 
not equations (1) and (2) are satisfied. 

-TJneway to check whether or not equations (1) and (2) have been met is t< 
between S a (t) and S b (t). Specifically one could compute the cross correlation 



f S a (t) S b (t) dt = \ { s a (t) s b (t) + n a (t) n b (t) ) dt + j ( s a (t) n b (t) + s b (t) n a (l) ) dt 

n a (t) and n b (t), may be 



which assuming that the signals, s a (t) and s b (t), are uncorrelated with the 
expressed as 



J S a (t) S b (t) dt = J ( s a (t) s b (t) + n a (t) n b (t) ) dt. 

A normalized cross correlation may be computed as follows 

I S a (t) S b (t) dt / ( s a (t) sb(t) + n a (t) n b (t) ) dt 



(J S2 a (t) dt I S2 b ( t ) dt]0-5 [J ( s 2 a(t) + „2 a{t) } dt f ( 5 2 b(t) + „2 b(t) } dt] 0.5' 



The solulions which give the extrema of this cross correlation functional are given by equation (2). Any 
input signals which are not expressible in the form given in equation (2) yield reduced normalized 
correlation values. In Figure 1, these extrema of the normalized correlation functional are plotted vs. the 

Copyright ©JMMasimo Corporation ™ 26052 Merit Circle Suite 103, Laguna Hills, CA, 92653 
^J"his document contains confidential and proprietary information. 



"Red" Signal - to - Noise ratio and W|/«>2. One should note that the normalized correlation is not 
necessarily unity and precaution must be taken, particularly when the difference between arterial 
saturation and the background saturation is large. 

One may rigorously prove that the magnitude of the normalized cross correlation is always less than or 
equal to unity by use of the inner product 

x.y =fx(t)y(0dt. 
Using the above notation, the normalized cross correlation may be written as 
V«b + Vb 



[(s a .s a + n a .n a )(s b .s b + n b .i Jb )]0-5 

or equivalent^ 

Sa-Sb 

[S a -s a s b .s 5] o- 5 ' 

But, according to the Schwarz Inequality 

(S a .S b )2 5 (S a .S a )(S b .S b ). 
Hence the normalized correlation is always between 1 and -1. 
ALGORITHM IMPLE MENTATION 
I. Compute the quantity 

£s a s b 



ES 2 a Zs2 b] 0-5" 
2. Compare quantity to a threshold, say 0.75 

THINGS THAT STILL NEED TO BE INVESTIGATED 

Investigate performance vs. snapshot size 
Investigate performance vs signal/noise correlation 
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Exhibit B 



j^'tfpC of**!' Janet Cerrefetionyil54B 
t neanjjFoHarlKks 



$5"Z1PEAK,C 
©PEAK.C 



2KB C Header fte 

3 KB c source fie 

4KB CSourcefie 

3 KB CHeaderfle 

4KB C Source fie 

■4KB CHeaderffc 

12 KB CSourcefie 

6 KB C Source Re 

2KB CHeaderfle 

3 KB CSourcefie 

3KB CHtaderfit 

to KB csourcefle 

3 KB CHeada-ffe 
5 KB CSourcefie 
4KB C Header fie 
12KB CSourcerte 

7KB CSourcefie 

2KB CHeaderfte 

4 KB csourcefie 
2KB CHeadB-fte 

4KB CSourcefie 

4KB csourcefie 

2 KB CHeederfte 
13KB CSourcefie 

SKB CSourcefie 

1KB CHeadcrrie 

4KB CHeadrfle 

2KB CHcodtrfle 

4KB csourcefie 

4KB CSourcefie 

7 KB Fit 

3KB csourcefie 

6KB CSourcefie 

3 KB CHeaderfle 
9 KB GOD He 

10 KB Chteaderne 

14KB csourcefie 

10 KB C Source fit 

0 KB DOT He 

0 KB Text Document 

6KB CSourcene 

30KB CSourcefie 



2/19/1993 11:27AM 
2/19/1993 11:40 AM 
2/13/1993 11:45 AH 
2/23/1933 12:53 PH 
2/23/1993 1:26 PM 
2/23/1993 3:1SPM 
2/23/1993 5:41PM 
2/23/1993 5:59 PH 
2/24/1993 3:02 PM 
2/24/19933:031*1 
2/25/1993 12:00 PH 
2/2S/1993 1:00 PM 
2/25/1993 4:26 PM 
2/23/1993 4:Z6PM 
2/26/1993 5 ;S0PM 
3/2/1993 10:09 AM 
3/2/1993 1:12 PM 
3/2/1993 2:45 PM 
3/2/1993 2:47 PH 
3/2/1993 2:47 PM 
3/2/1993 SiSliPH 
3/3/1993 S:03PN 
3/3/1993 5:03 PH 
3/3/1993 5:28 PH 
3/S/1993 2:21PM 
3/S/1993 2:54PM 
3/5/3993 3:06 PH 
3/5/1993 3:06 PM 
3/S/1993 3i27PH 
3/5/1993 3:37 PH 
3/6/1993 1:36 PH 
3/6/1993 1:41PM 
3/8/1993 B:53 AH 
3/6/1993 10:21 AM 
3/8/1993 12:47 PM 
3/B/1993 12:50 PM 
3/8/1993 l:04 PM 
3/B/1993 1:07 PM 
3/8/1993:1: 15 PH 
3/8/1993 2:06 PM 
3/8/1993 3:44 PH 
3/8/1993 3:45 PM 
3/B/1993 3:48 PM 
3/8/1993 3:49 PH 
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Fie and FoHer Taste 



L. 



3KB CHeadernc 
8KB c Header fie 
1KB CHeaderfle 

ikb CHeaderrie 

2KB CHeaderffe 

1KB CHeaderffe 
2KB CHeaderfte 
1KB CHeaderffe 
2KB c Header fie 
2KB C Heads-He 



. 2KB CHta derfle 

iPkaxare {Device- IndyandM^iM^ erfte 



©ZIVSNCC 
©Z1NANC.C 



7 KB CSOUTCBfle 

3 KB CSourcefle 

5 KB CSourcefle 

2KB C Header fie 

3 KB C Source ffe 

3 KB CHeaderffe 

3 KB CSourcefle 

8KB CSourcefle 

8KB CSourcefle 

1KB CHeaderffe 

2KB CHeaderffe 

6KB CSourcefle 

4KB CSourcefle 

23KB CSourcefle 

5KB CSourcefle 

1KB CSourcefle 

2KB CSourcefle 

1KB CSourcefle 

1KB CSourcefle 



3K3 
1KB 

1KB CSourcefle 
1KB CSourcefle 
1KB CSourcefle 
1KB CSourcefle 

2KB CSourcefle 
1 KB CSourcefle 
1KB CSourcefle 

1KB CSourcefle 

1 KB CSourcefle 



1/12/19922:43 PM 
7/21/1992 11:27 AM 
11/6/1992 5 :33 PM 
11/10/1992 11:29 AM 
12/8/19923:55 PM 
12/30/1992 4:53 PH 
12/30/1992 4:54 PM 
1/16/1993 2 :34 PM 
1/16/19932:38 PM 
1/18/19935:11 PM 
1/25/1993 2 :l+PM 
2/3/1993 4:29 PM 
2/4/1993 11 :22 AM 
2/9/1993 2I23PH 
2/9/1993 2:24 PM 
2/9/3993 2:28 PH 
2/12/19934:40 PM 
2/12/1933 4:52 PM 
2/16/19932:37 PM 
2/16/1993 2 ;38PM 
2/16/19932:45 PM 
2/16/1993 2:45 PM 
2/16/1993 2:50 PM 
2/17/1993 10:32AM 
2/17/1993 1 1:27 AM 
2/17/19932:30 PM 
2/18/19933:11 PM 
2/10/19935:45 PM 
2/18/19936:53 PM 
2/18/19936:56 PM 
2/18/1993 6:57 PM 
2/18/19936:58 PM 
2/18/1993 6:59 PM 
2/18/1993 7:02 PM 
2/18/1993 7 :04 PM 
2/18/1993 7 :06 PM 
2/18/1993 7 :08 PM 
2/18/1993 7 : 10 PM 
2/18A993 7:13PM 
2/18/1993 7: 17 PM 
2/18/19937:19 PM 
2/lS/1993 7;20PM 
2/1B/199 3 7:22 PM 
2/ 18/1993 7:26 PM 
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Copyright (c) Masimo Corporation (tm) 1993. 
File: vector. h 
Description: 

Collection of routines that work on vectors. Designed for 
Masimo Standard Library. 

Public Functions: 

min_index Returns index to minimum element in vector 

max_index Returns index to maximum element in vector 

innerjproduct Returns point by point product of vectors 

remove_mean Removes the mean value from each point in a vector 

Computes Dot product of two vectors 
Computes Correlation Coefficient 

Private Functions: 

Notes : 

History: 

01/15/93 
01/25/93 
02/25/93 
02/25/93 



RJM 
RJM 
RJM 
RJM 



Initial Version 
add remove_mean 
add VECTOR_DotProduct 
Add VECTOR CrossCoreO 



# include <masimo.h> 

extern INT32 

min_index ( 

FLOAT 3 2 *dat< 
INT32 size! 



extern INT32 
max_index ( 

FLOAT 3 2 

INT32 



*data, 
size) ; 



extern VOID 

inner_product ( 

FL0AT32 *vl, 

FL0AT32 *v2, 

INT32 size, 

FLOAT32 *v3) ; 

extern VOID 

remove_mean ( 

FLOAT32 *in, 

FLOAT32 *out, 

INT32 size) ; 

extern FLOAT 3 2 
VECTOR_DotProduct ( 

FLOAT32 *u, 

FLOAT32 *v, 

INT32 size) ; 

extern FLOAT3 2 
VECTOR_CrossCorr0 ( 

FLOAT 3 2 *u, 

FLOAT 3 2 *v, 

INT32 size) ; 



/* input, pointer to data 

/* input, number of data points 



/* input, pointer to data 

/* input, number of data points 



/* input, first vector 
/* input, second vector 
/* input, size of vectors 
/* output, v3 = vl * v2 



/* input, data 

/* output, data with mean removed 
/* input, number of points 



/* input, first vector 
/* input, second vector 
/* input, size of vectors 



/* input, first vector 
/* input, second vector 
/* input, size of vectors u and - 
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File: vector. c 

Description: 

Collection of routines that work on vectors. 
Masimo Standard Library. 

Public Functions: 
min_index 
max_index 
inner_jproduct 
remove_mean 
VECTOR_DotProduct 
VECTOR_CrossCorrO 



Returns index to smallest element in vector 
Returns index to largest element in vector 
Computes the inner product of two vectors 
Removes the mean value from each point in a vector 
Compute Dot Product of two vectors 
computes cross correlation time=0 



Private Functions: 



01/15/93 RJM Initial Version 
02/25/93 RJM Add VECTOR_DotProduct 
02/25/93 RJM Add VECTOR_CrossCorrO 



ftdefine MODULE_ID 

#include <masimo.h> 
#include <math.h> 
#include "vector. h" 



platform definitions 
sgrt function 
self 



INT32 

min_index ( 
FLO AT 3 2 
INT32 

INT32 
INT32 
FLOAT32 



*data, 
size) { 

i; 

j; 

tmp; 



input, pointer to data 
input, number of data points 



BUG1 (data) ; BUG1 (size) ; 
j = 0; 

tmp = data [j ] ; 

for (i = 1; i < size; i++) { 
if (data[i] < tmp) { 
tmp = data [i] ; 



} 



INT32 

max_index ( 

FLOAT32 *data, /* input, pointer to data */ 

INT32 size) { /* input, number of data points */ 

INT32 i ; 

INT32 j ; 

FL0AT32 tmp; 



BUG1 (data) ; BUG1 (size) ; 
j = 0; 

tmp = data [j ] ; 

for (i = 1; i < size; i++) { 
if (datafi] > tmp) { 
tmp = data [i] ; 

return ( j ) ; 



VOID 

inner jproduct ( 
FLOAT32 
PLOAT32 
INT32 
FLOAT32 

INT32 



*v2, 
size, 
*v3) { 



/* input, first vector 

/* input, second vector 

/* input, size of vectors 

/* output, v3 = vl * v2 



for (i = 0; i < size; i++) { 
V3 [i] = vl [i] * V2 [i] ; 

} 



VOID 

remove_mean ( 
FLOAT 3 2 
FLOAT 3 2 
INT32 



input , data 

output, data with mean removed 
input, number of points 



BUG1 (in) ; BUG1 (out) ; BUG1 (size) ; 
mean = 0.0; 

for (i = 0; i < size; i++) mean += in[i]; 
mean /= (FLOAT32) size; 

for (i = 0; i < size; i++) out [i] = in[i] 



/* sum all points */ 



FLOAT 3 2 

VECTOR_DotProduct ( 
FLOAT32 
FLOAT32 
INT32 



*u, 



size) { 



/* input, first vector 
/* input, second vector 
/* input, size of vectors 



BUG1 (u) ; BUG1 (v) ; BUG1 (size) ; 

sum = 0.0; 
130 for (i = 0; i < size; i++) { 

sum += u [i] * v [i] ; 

} 

return (sum) ; 



135 



} 



FL0AT32 
VECTOR_CrossCorrO ( 

FLOAT32 

FLOAT 3 2 

INT32 



*u, 



size) { 
den; 



/* input, first vector 
/* input, second vector 
/* input, size of vectors u and \ 



BUG1 (u) ; BUG1 (v) ; BUG1 (size) ; 



: SQRTF (VECTORJDotProduct (u, 
VECTOR_DotProduct (v, 



size) * 
size) ) ; 



if (den ==0.0) { 

c = 1.0; 
} else { 

c = VECTOR DotProduct (u, v, size) / den; 

} 

return (c) ; 
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Copyright (c) Masimo Corporation (tm) 1992, 1993. All Rights Reserved. 
File: stork. h 

Description: Stork header file included by ALL stork modules 
Public Functions: 
Private Functions: 
Notes : 



History: 




RJM 


11/06/92 


RJM 


12/08/92 


RJM 


12/16/92 


RJM 


01/04/93 


RJM 


01/05/93 


RJM 


01/14/93 


RJM 


01/15/93 


RJM 


01/16/93 


RJM 


01/18/93 


RJM 


01/18/93 


RJM 


01/19/93 


RJM 


01/25/93 


RJM 


01/25/93 



Initial version (vl.29) 
Prevent multiple inclusion of stork. h 
Change vsinc params to all floats 
Add MAX_BINS 

Add seed to bin structure 

Put conditional code around "C" only stuff 
Make memory pool sizes platform dependant 
Change sat result structure for filtering 
Make HR a vector 
Add more sat params 

Add hr_max_delta and sat_max_delta params 
Make SAT equation adjustable 
Make HR algorithms switchable 



#ifndef STORK 
#define STORK 

#define VERSION 1.541 /* floating point software version number */ 
#define VSINC_LEAN 0 /* Generate code for Lean Architecture */ 



#define ADC_SNAPSHOT_SIZE 



Memory Allocation 



#if (PLATFORM == 1) 

#define MEM_STACK_SIZE 
#define MEM_POOL_SIZE 
#define DEBUG_MEM_SIZE 
#define GDD_BUF_SIZE 
#define DBG_FIFO_SIZE 

#else 

#define MEM_STACK_SIZE 
#define MEM_P00L_SIZE 
ttdefine DEBUG_MEM_S I ZE 
#define GDD_BUF_SIZE 
ttdefine DBG FIFO SIZE 



32768 
512 

20000 
100 



48 
26000 
512 
6000 
100 



/* data memory 
/* data memory 
/* data memory 
/* pmda (wa 
/* pmda 



#endif 

/* The following contant controlls ALL debug code for the STORK 
/* beta software. No required code should be conditioned on this 
/* constant. 



#define DEBUG 1 



/* enables debug code 



*/ 



/' 



The following constants control what data is captured in the dbg_fifo. 
The dbg_fifo is controlled by the ADC interrupt routine. 



#define DEBUG_SIM 
#define DEBUG_ADC 
#define DEBUG_LED 
ttdefine DEBUG_SUMS 
ftdefine DEBUG_RAW 
ttdefine DEBUG_HSL 



/* simulate ADC data */ 

/* stores raw adc conversions in fifo */ 

/* stores tick mark for LED change */ 

/* stores the reading for each state */ 

/* stores raw ir and raw red */ 

/* stores High Speed Link data */ 



/ „ 

Low Speed Link Control 
*/ 

#define LSL_CHECK_OVERRUN 1 /* check HC11 FIFO */ 

#define LSL_GDD_GROUP_SIZE 6 /* number of gdd packets per lsl_tick */ 

/* 

High Speed Link Control 

*/ 

#define HSL_RECORD 0 /* records ADC data after ambient subtract */ 

#define HSL_PLAYBACK 0 /* uses data from HSL instead of ADC */ 

#define HSL_TEST_PATTERN 0 /* generate rotating l's pattern */ 

#define IMPULSE 0 /* generate impulses for usrl */ 

/* - - 

The following constants control the ANC. 

- — — - - */ 

#define ANC_FIXED_CELL 0 /* fixed cell normalized anc */ 

/* 

System Wide Type definitions 

*/ 

#define FIXED_CELL 1.0 /* anc_types */ 

#define AC CORD IAN 2.0 /* anc_types */ 

/* The following constants describe the size of the following */ 

/* structures. */ 

#define MAX_BINS 6 
#define NUM_HR_ALG 7 

#define BIN_RESULTS_SIZE 8 

idefine SAT_PARAMS_SIZE 22 

#define SAT_RESULTS_SIZE 14 

#define HR_PARAMS_SIZE 7 

#define HR_RESULTS_SIZE 8 

#ifndef ASSEMBLY 

typedef struct { 

INT32 num_jpoints; /* number of input points */ 

FLOAT32 anc_type; /* type of anc */ 

FLOAT 3 2 lambda; /* nANC parameter */ 

FLOAT 3 2 min_error; /* nANC parameter */ 

INT32 anc_settle_time; /* nANC parameter */ 





FLOAT3 2 


num cells; 


/ * 


nANC parameter 


j 

*/ \ 






peak limitl ; 


/* 


peak picker parameter 


*/ I 




FLOAT3 2 


peak epsilon; 


/* 


peak picker parameters 


*/ I 




FLOAT3 2 


*bp coef ; 


/* 


pointer to filter coef . 


*/ | 




INT32 


bp_size; 


/* 


number of filter coef. 


*/ 




FLOAT3 2 


search first ; 


/* 


first saturation value . 


*/ 




FLOAT3 2 


search step ; 


1* 


saturation step size. 


*/ I 






search size ; 


/* 


number of steps 


*/ f 






bin count ; 


/* 


number of bins 


*/ | 




„® „,„ 


cov lituit ; 


/* 


covariance limit 


*/ 1 






DVHS init ; 


/* 


initial value for DVHS 


*/ I 




FL0AT3 2 


DVHS0 init ; 


/* 


initial value for first DVHS 


*/ 1 




FL0AT3 2 




^* 


seconds 


*/ j 






DVAST delta; 


/* 


change for DVAST 


*/ j 




BLOAT 3 2 


DWST_tiraer ; 


/* 


seconds 


*/ 1 




FL0AT3 2 


DWST delta; 


/* 


change in DWST 


*/ I 




FLOAT32 


max_sat_delta ; 


/* 


max change in sat w/motion 


*/ i 




INT32 


sat_eqn; 


/* 


saturation equation number 


'', \ 


155 


FLOAT32 


dual_lambda ; 




dual ANCs as VSINC1 inputs 






FLOAT32 


correlation; 


/* 


bins OK only if >= corr. 


*/ § 




} SAT PARAMS; 








j 




typedef struct { 










FLOAT32 


ir_rms ; 


/* 


signal rms 


*/ i 




FLOAT32 


rd_rms ; 


/* 


signal rms 


*/ i 




FLOAT32 


cov; 


/* 


covariance 


*/ 




FLOAT32 


arterial__sat ; 


/* 


resulting arterial sat 


*/ \ 




FLOAT32 


venous_sat ; 


/* 


resulting venous sat 


*/ 


165 


FLO AT 3 2 


seed_sat ; 


/* 


seed saturation 


*/ I 




FLOAT 3 2 


numjpeaks ; 


/* 


peaks in power curve 


*/ i 




FLO AT 3 2 


correlation; 


/* 


raw dataset correlation 


*/ \ 




} BIN RESULTS 










170 


typedef struct { 










FLOAT 3 2 


Spa02_f inal; 


/* 


resulting arterial sat 


*/ ! 




FLOAT32 


Spv02_f inal; 


/* 


resulting venous sat 


* 1 \ 




FLOAT 3 2 


Spa02 unfil; 


/* 


unfiltered arterial sat 


Jl 1 




FLOAT32 


Spa02_fill; 


/* 


first stage of filtering 


*/ I 




FLOAT 3 2 


Spa02 dual; 


/* 


dual ANCs as input to VSINC1 


* 1 




FLOAT 3 2 


DVHS; 


/* 


default value high sat 


*/ I 




FLOAT32 


DVHS0 ; 


/* 


default value high sat first 


*/ | 




FLOAT32 


DVAST; 




default value art sat (PR) 


*/ \ 




FLOAT 3 2 


DVAST time; 


/* 


time DVAST was last updated 


*/ 1 
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FLOAT32 


DWST; 


/* 


default value ven sat (PR) 


*/ i 




FLOAT32 


DWST_time; 


/* 


time DWST was last updated 


*/ i 




FLOAT32 


art_omega_pr ; 


/* 


arterial omega for pulse rate */ | 




FLO AT 3 2 


ven_omega_pr ; 


/* 


venous omega for pulse rate 






FLOAT 3 2 


time; 


/* 


time for these results 


% 1 
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BIN RESULTS bins [MAX BINS] ; 


/* 


bin results 






} S AT_RE SULTS 








| 




typedef struct { 
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FLOAT 3 2 


resolution; 


/* 


HR resolution in BPM 


*/ j 




FLOAT 3 2 


raw_sampling_rate ; 


/* 


raw data sampling rate in Hz 






FLOAT32 


min_hr ; 


/* 


minimum hr allowed (BPM) 


*/ j 




FLOAT32 


max_hr ; 


/* 


maximum hr allowed (BPM) 


*/ ! 




FLOAT32 


nst; 


/* 


noise spec thresshold 


*/ s 
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FLOAT 3 2 


max_hr_delta; 


/* 


maximum change in HR w/motion*/ f 




INT32 


hr_enable; 


/* 


bits control with alg 






FLOAT 3 2 


hr_lambda ; 


/* 


ANC algo. uses this lambda 






} HR PARAMS; 
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typedef struct { 






1 




FLOAT 3 2 


hr [NUM HR ALG] ; 


/* 


(BPM) from each alg 


*/ j 




FLOAT32 


hr_final; 


/* 


(BPM) after all filtering 






FLOAT32 


BPM_jper_bin; 


/* 


BPM per bin of spectrum 






INT32 


num_bins ; 


/* 


number of bins per spectrum 


*/ 1 



205 HR_RESULTS; 
#endif 
#endif 



Exhibit G 



\ 



! 



Copyright (c) Masimo Corporation (tm) 1992, 1993. All Rights Reserved. 
File: stork. gdd 

Description: Debug Display Control File 
Public Functions: 
Private Functions: 



15 


Channels 97 to 88 


are for status line j 




History: 




i 




RJM 


12/10/92 


Initial version J 
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RJM 


12/21/92 


Add algorithm parameter | 




RJM 


12/22/92 


Add bin results table and plots j 




RJM 


12/28/92 


Add new HR algorithm results f 




RJM 


12/30/92 


Change to one algorithm | 




RJM 


12/30/92 


Add peak_epsilon parmaters | 


25 


RJM 


01/05/93 


Add Seed sat to bin data, remove others 




RJM 


01/16/93 


New SAT results structure J 




RJM 


01/18/93 


Add dvhs_init and dvhs0_init parameters 




RJM 


01/18/93 


Add hr5 result \ 




RJM 


01/18/93 


Add more sat params 7 | 


30 


RJM 


02/22/93 


Real-Time data Plotting ] 




RJM 


02/23/93 


Clean up Channel Assignments j 


35 


Debugl 


screen 





P 0,2 
DEBUG SCREEN 1 
30, 050, 025, 0375, 125, 7 
34,050,175,0375,275,5 
33,050,325,0375,425,5 
32,050,475,0375,575,5 
20,775,025,1000,100,5 
21,775,150,1000,225,0 
22,775,230,1000,305, 0 
23,775,310,1000,385, 0 
24,775,390,1000,465,0 
25, 775,470,1000,545, 0 
31,775,600,1000,725,5 
00,400,005,0000,000, 0 
01, 575, 005, 0000, 000, 0 
03,575,225,0000,000,0 
04,575, 370, 0000, 000, 0 
02,010,610,0000,000,0 
17,400,490,0014,000,12 
98, 000, 741, 0014, 000,12 
99, 000, 741, 0015, 000, 12 
97,000,755,0015,000,12 



display key = : 

raw data 
analysis data 
plethy 
freq spect 
power curve 
bin 1-5 



ADC Snapshot 
system params 
Sat Results 
HR Results 
Calib Results 
bin results 
ADC OOR count 
' Working ' Text 
Version# 
file name. . 



Debug2 screen 
P 1,1 

DEBUG SCREEN 2 
30,050,025,425,150,5 
32,050,225,425,330,7 
33,050,400,425,525,5 
00,450, 005, 000, 000, 0 
01, 625, 005, 000, 000, 0 



Raw, Spectrum, Plethy, tables 
; display key = 2 
; raw 
; plethy 

; System Parameters 
; Sat Results 



03,625,225,0000,000,0 
04,625,400,0000,000,0 
98,050,570,014,000,12 
99,050,570,015,000,12 
97,050,585,015,000,12 
END 



HR Results 
Calib Results 
■Working 1 Text 
Version* 
file name. . 



Debug3 screen power curve, bin graphs 



P 2,1 

DEBUG SCREEN 3 
20,050,020,380 
21,050,180,380 
22,050,330,380 
23,450,020,780 
24,450,180,780 
25,450,330,780 
02,060,460,000 
98, 050,590, 014 
99,050,590,015 
97,350,590, 015 



,130,7 
280,7 
430,7 

,130,7 
280, 7 
430, 7 
000,0 
000,12 
000,12 
000,12 



display key = 3 



bin results 
' Working 1 Text 
Version* 
file name. . 



Debug4 screen Raw, Spectrum, Plethy, tables 



P 3,0 
DEBUG 
30, 050, 025, 
32,050,180, 
33,050,330, 
00,300,-15, 
01,470,-15, 
03,470,195, 
04,470,320, 
98, 000,455, 
99, 000,455, 
97, 000,470, 



270,125,5 
270,280,5 
270,430,4 
000,000,0 
000,000,0 
0000,000,0 
0000,000,0 
014,000,12 
015,000,12 
015,000,12 



display key = 4 



f req 
plethy 

system params 
sat results 
HR Results 
Calib Results 
1 Working ' Text 
Version* 
file name. . 



Debug5 screen power < 



P 4, 0 
DEBUG , 
20,050 
21,050 
22,050 
23,375 
24, 375 
25, 375 
02, 020 
98, 020 
99, 020 
97,250 



SCREEN 5 

020,295,100,7 
,135,295,215,7 
250,295,320,7 
020,620,100,7 
135,620,215,7 
250,620,320,7 
345, 000, 000, 0 
,471, 014, 000, 12 
,471, 015, 000, 12 
471,015,000,12 



e,bin graphs 
display key = 



bin results 
1 Working ' Text 
Version* 
file name . . 



; Customer Screen #1 

P 5, 0 

Customer Screen 

52, 000,000,015,000,28 

77,050,050,450,150,21 

78, 050,200,450,300,21 

53, 500, 080, 014, 000, 16 

54, 500,230, 014, 000, 16 

15,500,110,014,000,16 

16,500,260,014,000,16 

31, 500,310, 630,400, 00 

17, 500,415, 011, 000, 08 



display key = 6 

stork legend 
Real time IR Data 
Real time RD Data 
Spa02 : 
Pulse : 
Spa02 value 
Pulse value 
ADC Snapshot 
; ADC OOR count HGK 



18,500,430,011,000,08 
14, 500,445, Oil, 000, 08 
40,500,460,012,000,08 



Correlation 
Tick 

calibrate: 1/0 message 



; Customer Screen #2 

P 6,0 

Customer Screen 

52,000,000,015,000,28 

35,050,050,450,150,21 

33,050,200,450,300,05 

53,500,080,014,000,16 

54,500,230,014,000,16 

15,500,110,014,000,16 

16,500,260,014,000,16 

31,500,310,630,400,00 

17,500,415,011,000,08 

18,500,430,011,000,08 

14,500,445,011,000,08 

40,500,460,012,000,08 

END 

Front End Screen 

P 7,1 

Front End Screen 

04, 500, 010, 0000, 000, 0 

31,500,200,775,575,7 

17,050,050,014,000,12 

77,050,200,450,375,23 

78,050,400,450,575,23 

40,050,020,011,000,28 



display key = 7 

stork legend 

raw analysis portion 

plethy 

Spa02 : 

Pulse : 

Spa02 value 

Pulse value 

ADC Snapshot 

ADC 00R count H6K 

Correlation 

Tick 

calibrate: 1/0 message 



Calib Results 
ADC Snapshot 
ADC OOR count 
Real time IR Data 
Real time RD Data 
calibrate: 1/0 message 



Channel Definitions 



Usrl snapshot 



G 30 

Raw Data 
Point number 
Amplitude 



G 31 

ADC SNAPSHOT 
Point Number 
Amplitude 



Frequency Spectrum 



G 32 

Spectrum 
HR (BPM) 
Amplitude 



Master Power Curve 



G 21 
Bin #1 
Sat 



G 22 
Bin #2 
Sat 
Amp 



G 23 
Bin #3 
Sat 
Amp 



Bin #4 

Sat 

Amp 



G 25 
Bin #5 
Sat 
Amp 

Plethysmograph 

G 33 

Clean Plethysmograph 

Time 

Amp 

; analysis data 

G 34 

Analysis Data 

Index 

Amp 

; raw analysis data 

G 35 

Raw Data (before Masimo) 

Index 

Amp 

Info Strings 

F 98,0 

Working 

F 99,0 

Stork Version %6.3f 

; Algorithm Parameters 

T 0,1,32 
ALG 1 

covJLimit, %9.1f 
num_j>oints , %9 . Of 



anc_type,%9.0f 
lambda, %9.4f 
min_error, %9 . 4f 
anc_settle,%9.0f 
num_cells, %9. Of 
peak_limitl, %9.2f 
peak_epsilon, %9 . 2f 
search_f irst, %9.2f 
search_step, %9 . If 
search_size, %9 . Of 
bin_count, %9.0f 
sampling rate, % 9. Of 
dvhs_init , %9 . If 
dvhs0_init,%9.1f 
DVAST_timer,%9.0f 
DVAST_delta,%9.1f 
DWST_timer , %9 . Of 
DWST_delta,%9.1f 
max_sat_delta, %9 . If 
max_hr_delta, %9 . If 
sat eqn,%9.0f 
hr enable, %9. Of 
Cal Target, % 9. 4 f 
Cal Err Band,%9.4f 
Cal Set. Time, %9. Of 
Cal IR PGA, %9.2f 
Cal RD PGA, %9.2f 
hr lambda, %9.4f 
dual lambda, % 9. 4 f 
Correlation, %9.4f 

Sat Results 

T 1,1,14 
ALG 1 

>Spa02_f inal , %9 . 4f 
>Spv02_f inal , %9 . 4f 
>Spa02_unfil, %9.4f 
>Spa02_fill,%9.4f 
>Spa02_dual,%9.4f 
DVHS, %9.4f 
DVHS0,%9.4f 
DVAST,%9.4f 
DVAST time,%9.4f 
DWST,%9.4f 
DWST time,%9.4f 
>art_omega_pr , %9 . 4f 
>ven_omega_pr , % 9 . 4 f 
>time, %9.4f 

HR Results 

T 3,1,8 
Alg 1 

>hr0 (ir),%9.4f 
>hrl (sig_ref ) , %9.4f 
>hr2 (deriv) ,%9.4f 
>hr3 (ir/n),%9.4f 
>hr4 (sig/n) ,%9.4f 
>hr5 (spec5) ,%9.4f 
>hr6 (ANC),%9.4f 
>hr_final,%9.4f 

Bin Results 

T 2,6,8 
Master Bin 
Bin 1 



Bin 2 
Bin 3 
Bin 4 
Bin 5 

ir_rms, %9 .4f 
rd_rms, %9 .4f 
>cov, %9.4f 
>arterial sat,%9.2f 
> venous sat , %9 . 2f 
>seed sat,%9.2f 
>num_p>eaks , %9 . Of 
>corr,%9.6f 

; Calib Results 

T 4,1,10 
C Results 
ir_trans,%9.4f 
rd_trans,%9.4f 
desired_ir,%9.4f tnA 
desired_rd, %9.4f raA 
actual_ir,%9.4f mA 
actual_rd,%9.4f mA 
ir_gain, %9.2f 
rd_gain,%9.2f 
ir_delta,%9.4f 
rd_delta, %9.4f 

; display strings 

M 40,0 
Calibrating 
S 52,1 

Masimo Stork 

S 53,0 

Spa02 : 

S 54,0 

Pulse : 

F 14,0 

Tick %12.0f 

F 15,0 

%3.0f %% 

F 16,0 

%3.0f BPM 

F 17,0 

ADC OOR Count %8.0f 
F 18,0 

Sensor Pos. %8.4f 

Real-Time Plotting 
R 77,1,400 

Plethysmographies Data (before Masimo} 
Time 

Amplitude 
R 78,1,400 

Plethysmographies Data (after Masimo) 
Time 

Amplitude 
X 12, 77, 0, 9 
IR 

X 13,78,0,12 



405 ; end of items 



Exhibit H 



1 /* 

Copyright (c) Masimo Corporation (tm) 1992, 1993. All Rights Reserved. 



File: iawa.c 

Description: 

Industry Algorithm Without Averaging 

Public Functions: 
init_iawa ( ) ; 
iawa_sat ( ) ; 

Private Functions: 

Notes : 



RJM 11/05/92 

RJM 11/20/92 

RJM 12/08/92 

RJM 01/05/93 



Based on s . c (VI. 29) 
use param and result structures 
Make compatible with VSINC 
Compute separate seed for each bin 



#define MODULE_ID 


1015 




#include 


<masimo.h> 


/* 


platform specification 


*/ 


# include 


<memory . h> 


/* 


special memory allocation 


*/ 


# include 


< debug. h> 


/* 


assert, grab 


*/ 


#include 


<math.h> 








# include 


" stork. h" 


/* 


system parameters 


*/ 


#include 


"fir.h" 


/* 


block fir module 


*/ 


#include 


"sateqn.h" 


/* 


saturation equations 


*/ 


#include 


"iawa.h" 


/* 


self 


*/ 


# include 


"vector .h" 


/* 


cross corr 


*/ 


#include 


"gddcli.h" 


/* 


General Debug Display Client 


#include 


" spt .h" 









*data, 
first, 
num, 
code) ,• 



#if ( PLATFORM == 1 
extern VOID 
presentation ( 

FLOAT 3 2 

INT32 

INT32 

INT32 

#endif 



VOID 

init_iawa (void) { 

/* nothing yet */ 



VOID 

compute_log ( 

FLOAT 3 2 *in, 
FLOAT32 *out, 
INT32 n) { 



/* this BS is to allow UI . C code to store. 



/* input data 
/* output data 
/* number of points 



/* This routine computes the natural log of the input signal */ 

/* Because of the filtering in the demodulator, there are times */ 

/* when the resulting signal can go below zero, this log routines */ 

/* limits small values so that the log is gracefully limited. */ 

#define LOG00001 -11.51293 



INT32 



BUG1 (in) ; BUG1 (out) ; BUG1 (n) ; 



for (i = 0; i < 
if (in[i] > 
else out [i] 

} 

75 



n; i++) { 

0.00001) out[i] = LOGF(in[i] ) ; 
= LOG00001; 



VOID 

iawa_sat ( 

FLOAT 3 2 *log_ir J /* input, normalized ir data */ 

FLOAT32 *log_rd, /* input, normalized rd data */ 

SAT_PARAMS *sp, /* input, parameters */ 

S AT_RE SULTS *sr) { /* output, results */ 

INT32 i; 

INT32 j ; 

1NT32 b; /* bin counter */ 

INT32 bin_size; /* point in a bin */ 

IEJT32 f_size; /* size of filtered data */ 

BIN_RESULTS *r; /* points to results for each bin */ 

FLOAT32 *f_ir; /* filtered ir */ 

FLOAT32 *f_rd; /* filtered rd */ 

FLOAT 3 2 *xaxis; /* for green marker lines */ 

BUG1 (log_ir) ; BUG1 (log_rd) ; BUG1 (sp) ; BUG1 (sr) ; 

assert (log_ir != NULL) ; 
assert (log_rd != NULL); 
assert (sp != NULL); 
assert (sr != NULL); 

f_size = sp->num_jpoints - sp->bp_size; 

bin_size = (f_size - sp->anc_settle_time) / sp->bin_count; 
assert (f_size > 0) ; 

mark_mem() ; 

f_ir = get_mem(f_size) ; 
f_rd = get_mem(f_size) ; 
xaxis = get_mem(1000) ; 

/* band-pass filter IR and RD signals */ 

blkfir (log_ir, sp->num_points, sp->bp_coef, sp->bp_size, f_ir) ; 

blkf ir (log_rd, sp->num_points, sp->bp_coef, sp->bp_size, f_rd) ; 

/* output analysis data */ 

b=0; /* build a dataset that demarks the bins divisions */ 

assert (f_size-sp->anc_settle_time < 1000); 
for(i=0; i < f size - sp->anc settle time; i++) 
{ 

xaxis [i]=(FLOAT32) (b-5000) ; 
if(!(i%bin size)) b=10000-b; 

} 

j = sp->anc_settle_time- 1; 
b = f_size - j ; 

GDDCLI_Chart (SPT_AnalysisData, 3 ,b, 0.0, 1.0); 
GDDCL I_Data (SPT_AnalysisData, b, &f_rd [ j ] ) ; 
GDDCLI_Data (SPT_AnalysisData,b, &f_ir [ j ] ) ; 
GDDCLl_Data (SPT_AnalysisData,b, xaxis) ; 

#if (PLATFORM == 1) 

presentation (f_rd, j , b, 0 ) ; 
#endif 



/* calculate correlations */ 

sr->bins[0] .correlation = VECTOR_CrossCorrO (&f_ir [ j ] ,&f_rd[j] , f_size-j) ; 
for (b = 1; b <= sp->bin_count; b++) 
r = & (sr->bins [b] ) ; 

r- correlation = VECTOR_CrossCorrO (&f_ir [j] , &f_rd [j ] ,bin_size) ,- 
j+= bin size; 

} 

/* Square both signals */ 

for (i = 0; i < f_size; i++) 

f_ir[i] *= f_ir[ij; 
f rd[i] *= f rd[i] ; 

} 

/* Integrate using binning */ 

sr->bins [0] . ir_rms = 0.0; 
sr->bins [0] .rd_rms = 0.0; 
j = sp->anc_settle_time - 1; 

for (b = 1; b <= sp->bin_count; b++) 

r = &(sr->bins [b] ) ; 



for (i = 0; i < bin size; i++) 
{ 

r->ir_rms += f_ir[j]; 
r->rd_rms += f_rd[j]; 

} 

/* update master bin */ 

sr->bins [0] . ir_rms += r->ir_rms; 
sr->bins [0] .rd_rms += r->rd_rms; 

r->ir_rms = SQRTF (r- >ir_rms / bin_size) ; 
r->rd_rms = SQRTF (r- >rd_rms / bin_size) ,- 

r->seed_sat = calc_sat (r->rd_rras / r->ir_rms, sp->sat_eqn) ; 

/* compute master bin results */ 
r = &{sr->bins [0] ) ; 

r->ir_rms = SQRTF (r->ir_rms / f_size) ; /* was bin_size HGK ' 

r->rd_rms = SQRTF (r->rd_rms / f_size) ; 

r->seed_sat = calc_sat (r->rd_rms / r->ir_rms, sp- >sat_eqn) ; 
f ree_mem ( ) ; 



Exhibit I 



1 /* 

Copyright (c) Masimo Corporation (tra) 1992. All Rights Reserved. 



File: satfil.c 

Description: Spa02 output filtering algorithms 

Public Functions : 
reset_satf il 
sat filter 



Resets all filters 
runs all sat filters 



Private Functions: 
max_diff_alg 



RJM 10/09/92 Initial version 

RJM 10/12/92 Allow multiple filtering streams 

RJM 11/06/92 Bring up to coding standards 

RJM 01/16/93 Merge binning algorithms, DVHS filters, etc 

RJM 01/18/93 If COV is low, and 1 peak don't report Sv02 

RJM 01/26/93 Limit final Spa02 by highest seed of all bins 

RJM 01/26/93 Compare max-diff values to Spa02_final 



#define MODULE_ID 

#include <masimo.h> 
#include 11 stork, h" 
#include "satfil.h" 
#include "sateqn.h" 
#include <memory.h> 
ttinclude <debug.h> 
ttinclude <vector.h> 
#include <math.h> 

VOID 

reset_satf il ( 

SAT_PARAMS *sp, 
SAT_RESULTS *sr) { 

BUGl(sp) ; BUGl(sr) ; 



/* platform description file */ 

/* system parameters */ 

/* self */ 

/* calc_omega */ 

/* Special memory routines */ 

/* assert, grab */ 

/* vector functions */ 

/* */ 



/* input, saturation parameters 
/* inout, saturation results 
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sr 


>DVHS = sp->DVHS init; 




sr 


>DVHS0 = sp->DVHS0_init; 




sr 


>Spa02_final = 0.0; 




sr 


>Spa02 unfil = 0.0; 




sr 


>Spa02 fill = 0.0; 


55 


sr 


>Spv02_final = 0.0; 




sr 


>Spa02 dual = 0.0; 




sr 


>DVAST = 0.0; 




sr 


>DWST = 0.0; 




sr 


>DVAST time = -sp->DVAST 


60 


sr 


>DWST_time = -sp->DWST_ 



This routine finds the two bins with the largest difference 
between arterial and venous sat and of these two returns 
the saturation which is closest to the previous sat. For 
more detail see Stork design note 27 . 



Select the bins that have a good arterial and venous 
saturation and order these bins based on the difference between 
the venous and the arterial saturation. 

Select the two bins with the greatest difference. If there 
are no bins resulting from step 1. 

Of these two bins select the one closest to the previous 
Spa02_final value. 

If there is no previous value, then return the largest value. 



95 



diff alg{ 










SAT PARAMS 


*sp, 


/* 


input, saturation parameters 


*/ 


SATJRESULTS 


*sr) { 


/* 


inout, saturation results 


*/ 


INT32 


i; 








BIN RESULTS 


*b; 


/* 


general bin ptr 


*/ 


INT32 


*good bins; 


/* 


array of good bin numbers 


*/ 


FLOAT32 


*diffs; 


/* 


array of differences 


*/ 


INT32 


num_good_bins ; 


/* 


number of good bins 


*/ 


FLOAT 3 2 


sat; 


/* 


return value 


*/ 


FLO AT 3 2 


satl ; 


/* 


arterial sat of bin with largest diff*/ 


FLOAT32 


sat2 ; 


/* 


sat of bin with 2nd largest diff. 


*/ 


BUG1 (sp) ; BUG1 (sr) ; 








assert (sp ! 


= NULL) ; 








assert (sr ! 


= NULL) ; 









mark_mem ( ) ; 

good_bins = get_mem (sp- >bin_count) ; 
diffs = get_mem(sp->bin_count) ,- 

/* select bins that have good arterial and venous sat */ 
/* R154 : bin must also have good correlation HGK */ 

num_good_bins = 0; 

for (i = 1; i <= sp->bin_count; i++) { 
b = &(sr->bins [i] ) ; 

if (b->arterial_sat > 0 && b->venous_sat > 0 && 
b- correlation >= sp- correlation) 

good_bins [num_good_bins] = i; 

dif f s [num_good_bins] = b- >arterial_sat - b->venous_sat; 
num_good_bins++ ; 



/* Find Two bins with largest difference */ 

switch (num_good_bins) { 

case 0: 

sat = sr->Spa02_unf in- 
break; 



case 1: 



i = max_index(dif f s, num_good_bins) ; 

sat = sr->bins [good_bins [i] ] .arterial_sat; 

break ; 

default : 

/* find bin with largest diff */ 

i = max_index (dif f s, num_good_bins) ; 

satl = sr->bins [good_bins [i] ] . arterial_sat ; 

/* find bin with second largest diff */ 

diffs[i] = 0.0; 

i = max_index(diffs, num_good_bins) ; 

sat2 = sr->bins [good_bins [i] ] .arterial_sat; 

/* If there is no previous Spa02 choose highest sat */ 

if (sr->Spa02_final ==0.0) { 
sat = satl; 

if (sat2 > satl) sat = sat2; 
} else { 

/* Otherwise, choose one closest to Spa02_final */ 
sat = satl; 

if (FABSF (satl - sr->Spa02_final) > 

FABSF (sat2 - sr->Spa02_final) ) sat = sat2; 

} 

break ; 



f ree_mem ( ) ; 
return (sat) ; 



sat_f ilter 

This routine takes the master power curve results and the 
five bins and computes an overall arterial and venous 
saturation based on the following algorithm. For more 
detail see Stork design note 27. 



VOID 

sat_f ilter ( 

SAT_PARAMS 
SATJRESULTS 



■Br) { 
*mb; 



BIN_RESULTS ' 
FLOAT 3 2 delta; 
INT32 i; 
FLOAT 3 2 max_seed; 



max sat for final filter output 
min sat for final filter output 



/* input, saturation parameters 
/* inout, saturation results 



/* master bin ptr 

/* maximum seed from all bins 



BUGl(sp) ; BUGl(sr) ; 



assert (sp != NULL); 



assert (sr != NULL) ; 

/* If master bin cov is less than limit then use its sat value */ 
/* otherwise, use the sat of the max_diff_alg of the bins. */ 

mb = &(sr->bins [0] ) ; 

if {mb->cov < sp->cov_limit) { 

sr->Spa02_unfil = mb->arterial_sat; 

/* Update DVHS */ 

if (mb->arterial_sat > sr->DVHS0) { 
sr->DVHS = mb->arterial_sat; 
sr->DVHS0 = sr->DVHS; 

} 

/* DVAST algorithm */ 

if (mb->arterial_sat > 0.0 && mb->arterial_sat < sr->DVHS) { 
sr->DVAST = mb->arterial_sat; 
sr->DVAST_time = sr->time; 

} 

} else { 

sr->Spa02_unfil = max_dif f_alg (sp, sr) ,- 
/* DWST algorithm */ 

if (mb->venous_sat > 0.0 && mb- >num_peaks > 1) { 
sr->DWST = mb->venous_sat; 
sr->DWST_time = sr->time; 

/* first stage of filtering (DVHS) */ 
sr->Spa02_fill = sr->Spa02_unfil; 

if (sr->DVHS < sr->Spa02_fill) sr->Spa02_f ill = sr->DVHS; 

/* second stage of filtering */ 
/* If there is motion and a previous value, limit the change */ 
/* to max_sat_delta . */ 

if (mb->cov > sp->cov_limit && sr->Spa02_f inal > 0.0) { 

delta = sr->Spa02_f ill - sr- >Spa02_f inal ; 

if (delta > sp- >max_sat_delta> delta = sp- >max_sat_delta; 

if (delta < -sp->max_sat_delta) delta = -sp->max_sat_delta; 

sr->Spa02_f inal += delta; 
} else { 

sr->Spa02 final = sr->Spa02_f ill; 

} 

/* Find largest seed from all bins */ 

max_seed = mb->seed_sat; 

for (i = 1; i <= sp->bin_count; i++) { 

if (sr->bins [i] .seed sat > max_seed) max_seed = sr->bins[i] . seed_sat 

} 

if (sr->Spa02_final < max_seed) sr->Spa02_f inal = max_seed; 

if (sr->Spa02_final > MAXJ3PA02) sr->Spa02_f inal = MAX_SPA02; 

if (sr->Spa02_final < MIN_SPA02) sr->Spa02_f inal = MIN_SPA02; 

/* Venous logic */ 

if (mb->cov < sp->cov_limit && mb->num_j>eaks < 2) { 
sr->Spv02_final = 0.0; 



} else { 

sr->Spv02 final = tnb->venous_sat; 

} 

/* Compute arterial and venous sat for use by Pulse Rate algorithms */ 
sr->art_omega_pr = 0.0; 

if (sr->time - sr- >DVAST_time < sp->DVAST_timer) { 

if (FABSF (sr->DVAST - sr- >Spa02_f ill) < sp->DVAST_delta) { 
sr->DVAST_time = sr->time; 

sr->art_omega_j?r = calc_omega (sr->DVAST, sp->sat_eqn) ; 
} else { 

if (sr->Spa02_fill > 0.0) { 

sr->art_omega_pr = calc_omega (sr->Spa02_f ill, sp->sat_egn) ; 



sr->ven_omega_pr = 0.0; 

if (sr->time - sr- >DWST_time < sp->DWST_timer) { 

if (FABSF (sr->DWST - mb- >venous_sat) < sp- >DWST_delta) { 
sr->DWST_time = sr->time; 

sr->ven_omega_pr = calc_omega (sr->DWST, sp->sat_eqn) ; 
} else { 

if (mb->venous_sat > 0.0 && 

tnb->cov > sp->cov_limit && mb- >num_peaks > 1) { 
sr->ven_omega_pr = calc_omega (mb->venous_sat, sp->sat_eqn) ; 

/* satfil */ 



310 #ifdef STANDALONE 

#include <stdlib.h> 
#include <stdio.h> 

315 VOID main (void) { 



} 

320 #endif 



Exhibit K 



1 /* 

Copyright (c) Masimo Corporation (tm) 1992, 1993. All Rights Reserved. 



This platform 

independant user interface provides the ability to start and stop 
the oximeter, change system parameters, and control the input data. 
There are minor differences in these features listed below: 



Feature 

Start/Stop 
Change system param 
LED calibration 
PAN through data 
Rerun same data 
Standard Mode 

Public Functions: 
init_ui 

user_interface 

Private Functions: 
store_raw 
storejparam 
store_re suits 
single_step 
test_pattern 



PC Platform 



Yes 
Yes 
Yes 



DSP Platform 

Yes 
Yes 
Yes 



Standard initialization 
Start the user interface 



send raw data via LSL 

send parameters via LSL 

send results via LSL 

Run a single snapshot of data 

send a test pattern via LSL 



Notes: 






History: 






RJM 


12/09/92 


Initial Version based on main.c 


RJM 


12/28/92 


New HR module vsi_hr 


RJM 


12/30/92 


New default values for parameters 


RJM 


12/30/92 


Add peak_epsilon parameter 


RJM 


01/18/93 


Add dvhs_init and dvhs0_init parameters 


RJM 


01/18/93 


Add hrfil module 


RJM 


01/18/93 


Keep track of time 


RJM 


01/25/93 


Set default for DVHS0 to 95 


RJM 


01/25/93 


Make SAT eqn parameter 


RJM 


01/25/93 


hr_enable param 


RJM 


02/02/93 


Add command for DSP Platform version 


RJM 


02/11/93 


Validate data after get_snapshot 


RJM 


02/12/93 


New default values for coefs 


RJM 


02/19/93 


Integrate file list feature from HK. 


RJM 


02/23/93 


Add Calibration parameters 


RJM 


02/24/93 


Improve ADC Snapshot 


HGK 


03/08/93 


move data_ptr code to usrl_c 



# define MODULE_ID 


1024 




# include 


<masimo.h> 


/* 


Platform description 


*/ 


# include 


<pckeybrd.h> 


/* 


key definitions 


*/ 


# include 


<memory . h> 


/* 


memory manager 


*/ 


iinclude 


<debug.h> 


/* 


virtual platform debugger 


*/ 


# include 


<math.h> 








# include 


"sateqn.h" 


/* 


saturation equations 


*/ 


# include 


" stork. h" 


/* 


System parameters 


*/ 


# include 


"adc .h" 


/* 


Demodulation module 


*/ 


# include 


"usrl .h" 


/* 


Filtering and Decimation 


*/ 


# include 


"calib.h" 


/* 


LED calibration 


*/ 


#include 


"vsinc.h" 


/* 


VSI noise cancellation module 


*/ 


ttinclude 


"satfil.h" 


/* 


Saturation Filters 


*/ 


#include 


"hrfil. h" 


/* 


Pulse Rate Filters 


*/ 





# include 


"lsl.h" 


/* 


low speed link module 


*/ 


70 


# include 


"iawa.h" 


/* 


Industry algorithm without averaging 


*/ 




# include 


"post .h" 


/* 


post band pass filter 


*/ 




# include 


"vsihr .h" 


/* 


VSI Heart Rate Module 


*/ 




# include 


"vsipl .h" 


/* 


VSINC Clean photoplethysmograph 


*/ 




#include 


11 gddcli .h" 






* / 


75 


# include 


» spt .h" 


/* 


Stork Packet Types 


*/ 




# include 


"vsidual .h" 


/* 


dual ANC sat version 


*/ 




# include 


"ui .h" 


/* 


Self 


*/ 




#if (PLATFORM == 1) 









#include <stdio.h> 
static BOOL presentation_f lag; 
#endif 



typedef struct 
FLOAT32 
FLOAT32 
FLOAT32 
FLOAT32 

} PARAM_DESC ; 



min_val ; 
max_val; 
inc_val ; 
def val; 



/* minimum value 

/* maximum value 

/* amount to change 

/* default value 



typedef struct { 



SAT PARAMS 


*vp; 


/* 


pointer to vsinc params 


SAT RESULTS 


*vr; 


/* 


pointer to vsinc results 


FLOAT 3 2 




/* 


pointer to power curve 


HR PARAMS 


*vhp ; 


/* 


pointer to Heart Rate params 


HR RESULTS 


*vhr; 


/* 


pointer to Heart Rate results 


FLOAT 3 2 


*f spec; 


/* 


frequency spectrum 


FLOAT 3 2 


*pleth; 


/* 


clean plethysmograph 


CALIB PARAMS 


*cp; 


/* 


pointer to calibration params 


CALIB RESULTS 


*cr; 


/* 


pointer to calibration results 



ALGORITHM_SPEC; 



The following constants control the user adjustable 
parameters . 



#define 
#define 
#def ine 
#define 
#define 
#def ine 
#define 
. #def ine 
#define 
#define 
#def ine 
#define 
#define 
ttdefine 
#def ine 
#define 
#define 
#define 
#define 
#define 
#define 
#define 
#def ine 
#define 
#define 



P_cov_limit 

P_num_points 

P_anc_type 

P_lambda 

P_min_error 

P_anc_settle_time 

P_num_cells 

P_jjeak_limitl 

P_jpeak_epsilon 

P_search_f irst 

P_search_step 

P_search_size 

P_bin_count 

P_sampling_rate 

P_dvhs_init 

P_dvhsO_init 

P_DVAST_timer 

P_DVAST_delta 

P_DWST_timer 

P_DWST_delta 

P_max_s a t_de 1 t a 

P_max_hr_de 1 1 a 

P_sat_egn 

P_hr_enable 

P_cp_target 



#define P_cp_error_bam 
#define P_cp_settle_tii 
#define P_cp_ir_gain 
#define P_cp_rd_gain 
#define P_hr_lambda 
#define P_dual_lambda 
#define P_correlation 



I_DESC pararas [N' 



/* 








min 


max 


inc 


ief */ 


/* 


cov_limit 


0 


*/ 


0.0, 


200.0, 


2.0, 


40 . 0 


/* 




1 


*/ 


100 .0, 


1000 . 0 , 


10 . 0 , 


650 . 0 


/* 




2 


*/ 


1.0, 


2.0, 


1.0, 


1 . 0 


/* 


Smbd?^ 


3 


*/ 


0.0, 


1.0, 


0 .000 


1, 0.75 


/* 


minerror 


4 


*/ 


0.0, 


1.0, 


0 .000 


1, O.OOC 


/* 




*/ 


0.0, 


1000.0, 


1.0, 


149 . 0 


/* 


num_cells 
peak_limitl 


6 


*/ 




50 . 0, 


1.0, 


12 . 0 


/* 


7 


*/ 


oo' 


1.0, 


0.01, 


0 . 02 


/* 


peak_epsilon 


8 


!/ 




1.0, 


0 .01, 


0 . 03 


/* 


search_f irst 


9 


*/ 


0.0, 


120.0, 


0.25, 


104 . 0 


/* 


search, step 




*/ 


0.0, 


5.0, 


0 .10, 


0 . 5 


/* 


search_size 


11 


*/ 


0.0, 


300 . 0, 


5.0, 


195 . 0 


/* 


bin_count 


12 


*/ 


1.0, 


5.0, 


1.0, 


5.0 


/* 


sampling rate 


13 


*/ 


10 . 0, 


200 . 0, 


1.0, 


62 . 5 


/* 


dvhs_init 


14 


*/ 


40.0, 


105.0, 


0.1, 


99 . 0 


/* 


dvhsO init 


15 


*/ 


40.0, 


105 . 0 , 


0.1, 




/* 


DVAST timer 


16 


*/ 


0.0, 


100 . 0 , 




10 ' 0 


/* 


DVAST delta 


17 


* / 






0l' 


2 . 0 


/* 


DWST timer 


18 


*/ 


0.0, 


100. o', 


i.o', 


10 .0 


/* 


DWST_delta 


19 


*/ 


0.0, 


100.0, 


0.1, 


4.0 


/* 


max sat delta 


20 


*/ 


0.0, 


20.0, 


0.1, 


2.0 


/* 


max_hr_delta 


21 


*/ 


0.0, 


40.0, 


0.1, 


4.0 


/* 


sat_eqn 


22 


*/ 


0.0, 


1.0, 


1.0, 


1.0 


/* 


hr_enable 


23 


*/ 


0.0, 




l:lk. 




/* 


cp_target 


24 


*/ 


0.0, 


i.o', 


0.25 


/* 




25 


*/ 


0.0, 


0.5, 


0.01, 


0.10 


/* 


cp_settle_time 


26 


*/ 


0.0, 


300.0, 


2.00, 


75.0 


/* 


cp_ir_gain 


27 




1:1: 


4096.0, 


0.1, 


2.0 


/* 


cp_rd_gain 


28 


*/ 






0.1, 


2.0 


/* 


hr_lambda 


29 






1.0, 


0.0001 


, 0.9999 


/* 


dual_lambda 


30 


*/ 


It 


1.0, 


0.0001 


, 0.9999 


/* 


correlation 


31 


*/ 


0.1, 


1.0, 


0.0001 


, 0.9999 



} 

/* ■ 



param_val [MAX_ALGORITHMS] [NUM_PARAMS] ; 



.t_ui (void) { 
/* nothing yet */ 



(PLATFORM == 1} 



presentation ( 

PLOAT32 *data, 
INT32 first, 



CHAR 
FILE 
INT32 



if ( !presentation_f lag) return; 
switch (code) { 

case 0 : printf ("Enter B/A filename { . B, .A will be added) 
scanf ( "%s" , filename) ; 
sprintf (str, "%s.B", filename) ; 
goto part2; 

sprintf (str, "%s.A" , filename) ; 

if ( (fp=f open (str, "w+" ) ) ==NULL) break; 
for (i=0;i<num;i++) { 

f printf (fp, "%f\n»,data[first+i] ) ; 

f close (fp) ; 
break; 



case 1 



} 

ttendif 



SAT P ARAMS 


*sp, 


/* 


input , 


parameters 


*/ 


FLOAT 3 2 


*ir, 


/* 


input , 


ir datav 


*/ 


FLOAT32 


*rd, 


/* 


input , 


red data 


*/ 


INT32 


size, 


/* 




number of points 


*/ 


INT32 


time_stamp) { 


/* 


input , 


starting point number 


*/ 



INT32 i,j,k; 

FLOAT32 *t; 

BUG1 (sp) ; BUG1 (ir) ; BUG1 (rd) ; BUG1 (size) ; BUG1 (time_stamp) ; 

#if (PLATFORM == 1) 

mark_mem ( ) ; 

t = get_mem(size) ; 

assert (t != NULL) ; 

j = sp->bp_size/2 + sp- >anc_settle_time; /* 1st pt */ 

k = size - sp->bp_size - sp->anc_settle_time; /* #usable pts */ 



( FL0AT3 2 ) t ime_s tamp , 1.0); 



for(i=0,-i<j ,-i++) t[i] = 
for( ;i<j+k ;i++) t [i] = 
for( ;i<size;i++) t [i] = 
GDDCLI_Chart (SPT_UsrlSnapShot , 3, size, 
GDDCLI_Data(SPT_UsrlSnapShot, size, ir) , 
GDDCLI_Data(SPT_UsrlSnapShot, size, rd) ; 
GDDCLI_Data(SPT_UsrlSnapShot, size, t ); 

GDDCLI_Chart (SPT_RawAnalysisData, 2,k, (FL0AT32 ) time_stamp, 1.0); 
GDDCLI_Data (SPT_RawAnalysisData, k, &ir[j]); 
GDDCLI_Data (SPT_RawAnalysisData, k, &rd[j]); 

free_mem() ; 
#endif 



static VOID 
265 store ^parameters ( 

INT32 num_algorithms) { 

INT32 i; 

270 BUG1 (num_algorithms) ; 

for (i = 0; i < num_algorithms; i++) { 



GDDCLI_Data (SPT_SystemParam, NUM_PARAMS , parara_val [i] ) ; 



} 

static VOID 
store_results ( 

ALGORITHM_SPEC 

INT32 

SATJPARAMS 
S AT_RE SULTS 
HR_PARAMS 
HR_RESULTS 
FLOAT32 
FLOAT32 
FLOAT32 
INT32 
INT32 
FLOAT32 



*alg, 

num_algorithms) { 



*vp; 
*vr; 
*vhp ; 
*vhr; 
*power ; 
*f spec; 
*pleth; 



/* Frequency SPECtrum 
/* clean pleth. 



temp ,- 

BUG1 (alg) ; BUG1 (num_algorithms) ,- 
mentry("store_results") ; 
/* frequency spectrum */ 
i = 0; 

vp = alg[i] .vp; 
vr = alg [i] .vr; 
power = alg [i] .power; 
vhp = alg[i] .vhp; 
vhr = alg[i] .vhr; 
fspec = alg [i] . f spec; 
pleth = alg[i] .pleth; 

/* frequency spectrum */ 

GDDCLI_Chart(SPT_FreqSpec, 1, vhr- >num_bins , vhp->min_hr, 

vhr- >BPM_j>er_bin) ; 
GDDCLI_Data (SPT_FreqSpec, vhr->num_bins, alg [i] .fspec) ; 

/* power curve */ 

GDDCLI_Chart (SPT_MasterPowerCurve, num_algorithms, vp->search_size, 

vp->search_f irst, - (vp- >search_step) ) ; 

for (i = 0; i < num_algorithms ; i++) { 

GDDCLI_Data (SPT_Master Power Curve, vp->search_size, alg[i] .power); 

for (j =0; j < vp->bin_count; j++) { 

GDDCLI_Chart ( SPT_PowerCurvel + j, num_algorithms , vp- >search_size, 

vp->search_first, - (vp->search_step) ) ; 

for (i = 0; i < num_algorithms; i++) { 

GDDCLI_Data (SPT_PowerCurvel + j, vp- >search_size, 

&(alg [i] .power [ ( j+1) *vp- >search_size] ) ) ; 



/* plethysmography signal */ 



GDDCLI_Chart (SPT_CleanPleth, 
for (i = 0; i < num_algorithms; 



algorithms, pleth size(vp), 0.0, 1.0); 

" H) { 



GDDCLI_Data(SPT_CleanPleth, pleth_size (vp) , alg [i] .pleth) ; 



#if (PLATFORM == 1) 

presentation (alg [0] .pleth, 0 , pleth_size (vp) ,1) ,- 
#endif 

/* store tabular results */ 

for (i = 0; i < num_algorithms ; i++) { 

GDDCLI_Data(SPT_SatResultS, SAT_RESULTS_SIZE, (FLOAT 3 2 * ) alg [i] . vr) ; 
^ GDDCLI_Data(SPT_HRResults, HR_RESDLTS_SIZE , (FLOAT32 *) alg [i] . vhr) ; 

for (i = 0; i <= alg [0] . vp- >bin_count ; i++) { 

GDDCLI_Data (SPT_BinResults , BIN_RESULTS_SIZE, 
^ (FLOAT32 *)&(alg[0] .vr->bins [i] ) ) ; 

/* master bins' correlation sent separately (times 100) */ 

temp = alg [0] .vr->bins [0] .correlation * 100.0; 
GDDCLI_Data (SPT_correlation, 1, &temp) ; 

mexit ("store results"); 

} 

static INT32 
compute_snapshot_size ( 

ALGORITHM_SPEC *alg, 

INT32 num_algorithms) { 

INT32 i; 
INT32 size; 

BUG1 (alg) ; BUG1 (num_algorithms) ; 

/* find largest snapsnot size */ 
size = 0; 

for (i = 0; i < num_algorithms; i++) { 

if (alg [i] . vp->num_points > size) size = alg [i] .vp->num_jpoints; 

return (size) ; 

} 

static VOID 
newjpatient ( 

ALGOR I THM_S PE C *alg, 

INT32 num_algorithms) { 

INT32 i; 

BUG1 (alg) ; BUG1 (num_algorithms) ; 

for (i = 0; i < num_algorithms; i++) { 
reset_satf il (alg [i] .vp, algti] .vr) ; 
reset_hrf il (alg [i] .vhp, alg [i] .vhr) ; 



static VOID 
single_step( 

ALGORITHM_SPEC *alg, /* input, algorithm structure 

INT32 num_algorithms, 

INT32 raw_snapshot_size) { 



INT32 
FLOAT32 



/* algorithm 
/* raw ir data 



*/ 
*/ 





FLOAT 3 2 


*rd; 


/* raw rd data 


*/ 


410 


FLOAT32 


* log_ir ; 


/* log of ir data 


*/ 




FLOAT32 


*log_rd; 


/* log of rd data 


*/ 




FLOAT32 


*adc_snap; 


/* snapshot of adc data 


*/ 




SAT PARAMS 


*vp; 






415 


SAT RESULTS 


*vr; 








HR PARAMS 


*vhp; 








HR RESULTS 


*vhr ; 








FLOAT 3 2 


*power ; 








FLOAT32 


*f spec; 


/* Frequency SPECtrum 


*/ 


420 


FLOAT32 


*pleth; 


/* clean pleth. 


*/ 




INT32 


time; 








BUGl(alg) ; 


BUG1 (num_algorithms) ; 


BUG1 <raw_snapshot_size) ; 




425 


mentry ( " single step 11 ) ; 








mark_metn ( ) ; 










/* Allocate data required for current algorithms */ 





ir = get_mem 

rd = get_mem 

log_ir = get_mem 

log_rd = get_mem 

adc_snap = get_mem 



( raw_snapshot_s i ze ) 
( r aw_snapshot_s i ze ) 
( raw_snapshot_si ze ) 
( raw_snapshot_size ) 
( ADC_SNAPSHOT_SIZE ) 



ADC_GetSnapshot (adc_snap) ; 
adc_snap [0] = -1.0; 
adc_snap[l] = 1.0; 

GDDCLI_Chart(SPT_ADCSnapShot, 1, ADC_SNAPSHOT_SIZE, 0.0, 1.0); 
GDDCLI_Data(SPT_ADCSnapShot, ADC_SNAPSH0T_SIZE, adc_snap) ; 

GDDCLI_SendPacket (SPT_ADC_00R, ( FLOAT3 2 ) ADC_Get0ut0f RangeCount ( ) ) 

time = get_snapshot ( ir , rd, raw_snapshot_size) ; 

/* validate this snapshot of data */ 

for (i = 0; i < raw_snapshot_size ; i++) { 

if (ir[i] > 2.0) assert (0); 

if (ir[i] < -2.0) assert(O) ; 

if (rd[i] > 2.0) assert (0); 

if (rd[i] <-2.0) assert(O); 

} 

store_raw(alg[0] .vp,ir, rd, raw_snapshot_size, time) ,- 

compute_log(ir, log_ir, raw_snapshot_size) ; 
compute_log (rd, log_rd, raw_snapshot_size) ; 

for (i = 0; i < num_algorithms ; i++) { 

/* setup pointers for current algorithm */ 

vp = alg [i] .vp; 
vr = alg [i] .vr; 
vhp = alg [i] .vhp; 
vhr = alg[i] .vhr; 

/* allocate variable sized results arrays */ 

alg [i] .power = get_mem (vp->search_size * (vp->bin_count + 1) ) 
alg [i].f spec = get_mem (vsi_spectral_size (vhp) ) ; 
alg[i] .pleth = get_mem (plethjsize (vp) ) ; 

power = algti] .power; 



fspec = alg [i] . fspec; 
pleth = alg [i] .pleth; 

/* time stamp data */ 

vr->time = (FLOAT32) time / vhp->raw_sampling_rate; 

iawa_sat ( log_ir , log_rd, vp, vr) ,- 

vsinc(ir, log_ir, log_rd, vp, vr, power); 
sat_f ilter (vp, vr) ,- 

vsi_hr( log_ir, log_rd, vp, vr, vhp, vhr, fspec); */ 

vsi_hr(ir,log_ir, log_rd, vp, vr, vhp, vhr, fspec); 

hr_f ilter (vhp, vhr, vp, vr) ; 

vsinc_pl(ir, log_ir, log_rd, vr->art_omega_jpr, vp, pleth); */ 
vsinc_jpl (ir, log_ir, log_rd, 

vr- >art_omega_pr , 
vhp->hr_lambda, /* HR's lambda for plethy. */ 

vp, pleth) ; 

vr- >Spa02_dual = 0.0; 

VSlDUAL_sat (vr->art_omega_jpr, 
vp- >dual_lambda , 
vp, 

ir, rd, log_ir, log_rd, 
&vr->Spa02_dual) ; 

GDDCLI_SendPacket(SPT_time_tick, (FLOAT32) time) ; 
GDDCLI__SendPacket (SPT_sat_f inal , vr- >Spa02_f inal) ; 
GDDCLI_SendPacket (SPT_hr_f inal, vhr- >hr_f inal) ; 

} /* for */ 

store_jparameters (num_algorithms) ; 
store_results (alg, num_algorithms) ; 

GDDCLI_NewLine ( ) ; 

GDDCLI_FlushBuffer() ; /* send data */ 

f ree_mem ( ) ; 

mexit ( "single_step" ) ; 



static VOID 
test_pattern (void) { 

/* send a test pattern to verify graphing stuff */ 



static VOID 

set_alg ( 

ALGORITHM_SPEC *alg, 

INT32 num_alg) { 

SAT_PARAMS *vp; 
HR_PARAMS *vhp; 
CAL IB_PARAMS * cp ; 

INT32 i; 

BUG1 (alg) ; BUG1 (num_alg) ; 

for (i = 0; i < num_alg; i++) 



vp = alg [i] . vp; 

vp- >cov_limit 

vp - >num_point s 

vp->anc_type 

vp->lambda 

vp- >min_error 

vp- >anc_settle_time 

vp- >num_cells 

vp- >peak_limitl 

vp- >peak_epsilon 

vp->bp_coef 

vp->bp_size 

vp- >search_f irst 

vp- >search_step 

vp->search_size 

vp- >bin_count 

vp->DVHS_init 

vp- >DVHS0_init 

vp- >DVAST_timer 

vp->DVAST_delta 

vp - >DWST_timer 

vp->DWST_delta 

vp- >max_sat_delta 

vp->sat_egn 

vhp = alg[i] .vhp; 

vhp- >resolution 

vhp- >raw_sampling_rate 

vhp->min_hr 

vhp->max_hr 

vhp->nst 

vhp- >max_hr_delta 
vhp - >hr_enabl e 

vhp- >hr_lambda 
vp- >dual_lambda 
vp- >correlation 

cp = alg [i] . cp; 
cp- >rd_min_mA 
cp- >rd_max_mA 
cp->ir_min_mA 
cp->ir_max_mA 
cp->target 
cp- >error_band 
cp- >ir_drive_l 
cp->rd_drive_l 
cp - > i r_dr i ve_2 
cp- >rd_drive_2 
cp- >settle_time 
cp->ir_gain 
cp- >rd_gain 
cp- >mu 



= param_val [i] 
= { INT3 2 ) param_ 
= (INT32)param_ 
= param_val [i] 
= param_val [i] 
= (INT32)param 
= (INT32) param_ 
■■ param_val [i] 
= param_val [i] 
= post_coef; 
: post_size; 
= parara_val [i] 
= param_val [i] 
= (INT32) param 
= (INT32 ) param 
■■ param_val [i] 
= param_val [i] 
■. param_val [i] 
= param_val [i] 
= param_val [i] 
= param_val [i] 
= param_val [i] 
= (INT32)param_ 



[P_cov_lirait] ; 
i_val [i] [P_num_points] ; 
_val [i] [P_anc_type] ; 
[P_lambda] ; 
[P_min_error] ; 
val [i] [P_anc_settle_time] ; 
_val [i] [P_num_cells] ; 
[P_peak_limitl] ; 
[P_jpeak_epsilon] ; 



[P_search_f irst] ; 
[P_search_step] ; 
_val [i] [P_search_size] ; 
_val [i] [P_bin_count] ; 
[P_dvhs_init] ; 
[P_dvhsO_init] ; 
[P_DVAST_timer] ; 
[P_DVAST_delta] ; 
[P_DWST_timer] ; 
[P_DWST_delta] ; 
[P_max_sat_delta] ; 
yal [i] [P_sat_egn] ; 



■■ 1.0; 

■ (INT32)param_val [i] [P_sampling_rate] ; 
■■ 25.0; /* bpm */ 

. 225.0; /* bpm */ 

■ 0.04; /* thresshold */ 
■■ param_val [i] [P_max_hr_delta] ; 

■■ (INT32)param_val [i] [P_hr_enable] ; 

■ param_val [i] [P_hr_lambda] ; 

: param_val [i] [P_dual_lambda] ; 

■ param_val [i] [P_correlation] ; 



. 20.0; /* mA */ 

■ 50.0; /* mA */ 
: 5.0; /* mA */ 

■ 50.0; /* mA */ 

■ param_val [i] [P_cp_target] ; 

■ param_val [i] [P_cp_error_band] ; 

5.0; /* mA */ 

■ 20.0; /* mA */ 
: 10.0; /* mA */ 

25.0; /* mA */ 

■ (INT32)param_val [i] [P_cp_settle_time] ; 
: param_val [i] [P_cp_ir_gain] ; 

■■ parara_val [i] [P_cp_rd_gain] ; 
0.1; 



#define DATA_STEP 57 /* amount to scroll : 
VOID 

user_interf ace ( 

INT32 num_algorithms) { 

ALGORITHM SPEC *alg; 



INT32 alg_num; /* current algorithm */ 





INT32 


param_num; 


/* 


current parameter number 


*/ 


615 


BOOL 
BOOL 
INT32 


quit; 

automation; 
lsl_cmd; 


/* 
/* 
/* 


quit program 

quit automated session? 

Low Speed Link command 


*/ 
*/ 
*/ 


620 


FLOAT32 
FLOAT32 
INT32 


*ir; 
*rd; 

preview_size 


/* 
= 650 


for data preview 


*/ 




INT32 


ss_size; 


/* 


snapshot size 


*/ 


625 


INT32 
INT32 
BOOL 
INT32 
FLOAT32 


i; 

j; 

cal_ok; 
time ; 
f ; 









BUG1 (num_algorithms) ; 
assert (num_algorithms • 



= MAX_ALGORITHMS) ; 



markjnem ( ) ; 

alg = get_mem{num_algorithms * sizeof (ALGORITHM_SPEC) ) ; 

/* now allocate structures for each algorithm */ 



i < num_algorithms ; i++) { 
vp = get_mem ( sizeof (SAT_PARAMS) ) ; 
vr = get_mem (sizeof (SAT_RESULTS) ) ; 
vhp = get_mem (sizeof (HR_PARAMS) ) ; 
vhr = get_mem (sizeof (HR_RESULTS) ) ; 
cp = get_mem( sizeof (CALIB_P ARAMS) ) ; 
cr = get_mem ( sizeof (CALIB_RESULTS) ) ; 



for (i = 0 
alg[i] 
alg[i] 
algti] 
alg[i] 
algti] 
algti] 



/* init param structures */ 

for (i = 0; i < num_algorithms; i++) { 
for (j = 0; j < NUM_PARAMS; j++) { 

param_val [i] [j] = paramstj] .def_val; 



param_num = 0; 
alg_num = 0; 
quit = FALSE; 

ir = get_mem (preview_size) ; 
rd = get_mem(preview_size) ; 



set_alg(alg, num_algorithms) ; 
newjpatient (alg, num_algorithms) ; 

store_parameters (num_algorithms) ; /* display initial parameters */ 
GDDCLI_FlushBuffer() ; 

CALIB_Init(alg[0] .cp, alg[0] .cr) ; 
do { 

GDDCLI_SendPacket (SPT_CalibrateMsg, 1.0); 
GDDCLI_FlushBuffer() ; 

cal_ok = CALIB_TwoPointCal(alg[0] .cp, alg[0].cr); 
GDDCLI_Data (SPT_CalibResults , CALIB_RESULTS_SIZE, 

(FLOAT32 *) (alg[0] .cr) ) 

} while (cal_ok == FALSE) ; 
GDDCLI_SendPacket (SPT_CalibrateMsg, 0.0); 



6DDCLI_FlushBuf f er ( ) ; 

time = get_snapshot (ir, rd, preview_size) ; 

store_raw(alg [0] .vp, ir, rd, preview_size, time); 
685 GDDCLI_FlushBuf f er ( ) ; 

presentation_f lag = FALSE; 

while (quit == FALSE) { 

690 

GDDCLI_HighLight (SPT_SystemParatrt, param_num, alg_num, 

param_val [alg_num] [param_num] ) ; 

GDDCLI_FlushBuf f er ( ) ; 
lsl_cmd = lsl_read(); 
695 GDDCLIJNormal (SPT_SystemParam, param_num, alg_num, 

param_val [alg_num] [param_num] ) ; 

GDDCLI_FlushBuf f er { ) ; 

switch (lsl_cmd) { 
700 case NLT: /* left arrow */ 

if (alg_num > 0) { 
alg_num -= 1; 

break; 

705 

case NRT: /* right arrow */ 

if (alg_num < num_algorithms-l) { 
alg_num += 1; 

710 break; 

case NUP: /* up arrow */ 
if (param_num > 0) { 
param_num -= 1; 

715 } 

break; 

case NDN: /* down arrow */ 

if (param_num < NUM_PARAMS-1) { 
720 param nura += 1; 

} 

break; 
case NINS : 

725 f = params [param_num] . inc_val * 100.0; /* HGK */ 

goto kplus; 

case NPLUS: 

f = params [param_num] . inc_val ; /* HGK */ 

73 0 #if (PLATFORM ==1) 

if (shift_key() ) f*= 10.0; 

#endif 
kplus : 

735 param__val [alg_num] [param_num] += f; 

if (param_val [alg_num] [param_num] > params [param_num] . max_val ) 
param_val [alg_num] [param_num] = params [param_num] . max_val ; 
break ; 

740 case NDEL : 

f = params [param_num] . inc_val * 100.0; /* HGK */ 

goto kminus; 

case NMINUS: 

745 f = params [param_num] . inc_val ; /* HGK */ 

#if (PLATFORM == 1) 

if (shift_key() ) f*= 10.0; 

#endif 



param_val [alg_num] [param_num] -= f; 

if (param_val [alg_num] [param_num] < params [paramjaum] . min_val ) 
param_val [alg_num] [param_num] = params [param_num] .min_val; 
break; 



case NRTUP: /* page up */ 

move_snapshot (-DATA_STEP) ; 

time = get_snapshot (ir, rd, preview_size) ; 

store_raw(alg [0] .vp,ir, rd, preview_size, time); 

GDDCLI_FlushBuffer() ; 

break; 

case NRTDN: /* page down */ 

move_snapshot { DATA_STEP) ; 

time = get_snapshot (ir, rd, preview_size) ; 

store_raw(alg [0] .vp, ir, rd, preview_size, time) ; 

GDDCLIJFlushBuffer () ; 

break ; 

case SCR: /* continuous operation */ 
set_alg(alg, num_algorithms) ; 

ss_size = compute_snapshot_size (alg, num_algorithms) ; 
/* calibrate2 (system_params.cal_target) ; */ 
reset_f ile_num() ; 
new_file(0) ; 

automation = TRUE; /* 'continue automated session' */ 

do 

{ 

time=0; 

new_patient (alg, num_algorithms) ; 
do 

{ 

single_step (alg, num_algorithms , ss_size) ; 
move_snapshot (USRl_AUTO_STEP) ; 

/* mating get_snapshot ( ) is in single_step. . */ 
if (lsl_status() ) 

automation = FALSE; /* discontinue session */ 

while (automation && ! end_of _data ( ) ) ; 
if (automation) 

GDDCLI NewFile ( ) ; 

} 

while (automation && ! end_of_f iles () ) ; 
break; 

case SCS: /* single step */ 

set_alg(alg, num_algorithms) ,- 

ss_size = compute_snapshot_size (alg, num_algorithms) ; 

single_step (alg, num_a 1 go ri thins , ss_size) ; 

break; 

case SCC: /* "C" calibrate command */ 
break; 

case SCT: /* "T" test pattern command */ 
test_pattern() ; 
break; 

case NESC: /* escape */ 
quit = TRUE; 
break ; 

case SCF: /* "F" new file command */ 
new_file(l) ; 



time = get_snapshot (ir, rd, preview_size) ; 

store_raw(alg [0] .vp,ir, rd, preview_size / time); 
GDDCLI_FlushBuf f er ( ) ; ; 
break; 

te SCB: /* "B" back file */ 
new_file(-l) ; 

time = get_snapshot (ir, rd, preview_size) ; 

store_raw(alg [0] .vp,ir, rd, preview_size, time) ; 
GDDCLI_FlushBuffer() ; 
break ; 



param_val [alg_num] [P_num_cells ]=2; 
param_val [alg_num] [P_search_size] =25; 
break ; 

#if (PLATFORM == 1) 

case SCW : presentation_f lag = 1- presentation_f lag; 
break; 

#endif 

default: /* ignore unknown commands */ 



} /* switch */ 
} /* while */ 



f ree_mem ( ) ; 
} /* user_interf ace : 



